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Abstract 

We present a way to deal with dispersion-dominated "shock-type" transition in the 
absence of completely integrable structure for the systems that one may character- 
ize as strictly hyperbolic regularized by a small amount of dispersion. The analysis 
is performed by assuming that, the dispersive shock transition between two different 
constant states can be modelled by an expansion fan solution of the associated modula- 
tion (Whitham) system for the short-wavelength nonlinear oscillations in the transition 
region (the so-called Gurevich - Pitaevskii problem). We consider as single- wave so 
bi-directional systems. The main mathematical assumption is that of hyperbolicity 
of the Whitham system for the solutions of our interest. By using general proper- 
ties of the Whitham averaging for a certain class of nonlinear dispersive systems and 
specific features of the Cauchy data prescription on characteristics we derive a set of 
transition conditions for the dispersive shock, actually bypassing full integration of 
the modulation equations. Along with model KdV and mKdV examples, we consider 
a non-integrable system describing fully nonlinear ion-acoustic waves in collisionless 
plasma. In all cases our transition conditions are in complete agreement with previous 
analytical and numerical results. 
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1 Introduction 



It is well known that the resolution of breaking singularities in dispersive media occurs 
through generation of short-wavelength nonlinear oscillations. The wave-like transition be- 
tween two smooth or constant hydrodynamic states is generally called a dispersive shock (or 
an undular bore, especially in the context of water waves). The main observable feature of 
the dispersive shock is formation of solitary waves in the vicinity of one of its edges. At 
the opposite edge, the wave structure degenerates into linear wavepacket. The phenomenon 
of the dispersive shock formation is quite ubiquitous and its physical contexts range from 
gravity water waves and bubbly fluid dynamics to space plasma physics, fibre optics and 
Bose-Einstein condensates. 

One should be clear from the very beginning that purely dispersive resolution of a shock is 
a physical idealization and some amount of dissipation is inevitably present in real continuous 
media so we distinguish between two types of dispersive shocks depending on the actual 
role of dissipation relative to that of nonlinearity and dispersion in the wave development 
and evolution. There is some terminological confusion in the literature where the terms 
"dispersive shock" , "undular bore" and "collisionless shock" are used as for the dispersion- 
dominated waves so for the waves where dispersion and dissipation are in balance. In both 
cases, dispersion plays the decisive role in the formation of local oscillatory structure but 
the global properties of the oscillations zone in the two cases are drastically different. 

The weakly dissipative undular bores (we use here this term to distinguish from our 
main subject - conservative dispersive shocks) despite their oscillatory structure exhibit 
global properties characteristic for classical, turbulent bores or shock waves: they have 
steady (though oscillatory in space) profile and constant effective width such that the speed 
of the shock propagation and the transition conditions can be derived within the frame of 
the classical theory of hyperbolic conservation laws (see for instance). The qualitative 
theory of such steady undular bores has been first developed by Benjamin and Lighthill [2] 
in the context of shallow-water waves and by Sagdeev j3j for rarefied plasma flows. The 
quantitative description of the weakly dissipative undular bores has been made in 
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jE] on the basis of the unidirectional KdV-Burgers equation and in |7] using the integrable 
version of the bi-directional Boussinesq equations modified by small viscous term. 

The developed in [2] [7] theory of steady undular bores, however, is valid only for the 
fully established regime when nonlinearity, dispersion and dissipation are in balance. Con- 
trastingly, in the dispersion- dominated case, the traditional analysis of the mass, momentum 
and energy balance across the undular bore transition can not be applied (at least directly). 
The reason for that is that the boundaries of the dissipationless undular bore (a dispersive 
shock) diverge with time, i.e instead of a single shock speed defined by the jump condi- 
tions one has now two different speeds s + > s~ determining motion of the transition region 
boundaries. These speeds , however, can not be found without the analysis of the nonlinear 
oscillatory structure of the transition region since in the dissipationless case dispersion not 
only dramatically modifies the fine structure of the shock transition but also, along with 
nonlinearity, determines its location. The "conservative" dynamics of dispersive shocks is of 
a considerable interest on its own and also, in many cases, can be viewed as an unsteady in- 
termediate asymptotic in a general setting when the small dissipation is taken into account. 
Physical examples of such expanding dispersive shocks include atmospheric undular bores 
(morning glory) |5], optical shocks in the long-distance optical communication systems [§| 
and collisionless shocks in Bose-Einstein condensates [TU] . 

In a weakly nonlinear case, when the original conservative system can be approximated by 
one of the exactly integrable equations, the study of the dispersive shock phenomenon have 
stimulated discovery of a whole new class of mathematical problems which can be broadly 
described as the singular semi-classical limits in integrable systems. In the Lax-Levermore- 
Venakides theory [TT) . [T2] (see also [T3j) developed originally for the KdV equation, and 
more recently extended to the defocusing NLS [Hj and focusing mKdV [THj equations the 
evolution of the dispersive shock is modelled by the zero dispersion limit of the exact multi- 
soliton (multi-gap) solution of the original dispersive wave equation. The main characteristic 
feature of the zero- dispersion limit is co-existence of smooth and rapidly oscillating regions 
in the solution after the breaktime. In the smooth regions the zero-dispersion limit exists in 
a strong sense and is given by the classical solution of the dispersionless equation whereas 
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in the oscillating regions this limit exists in a weak, averaged sense and is governed by a 
certain system of quasi-linear equations different from the dispersionless limit. This system 
turned out to coincide with the modulation equations obtained in 1960's by Whitham ^H] by 
averaging nonlinear single-phase wavepackets and later generalised by Flaschka, Forest and 
McLaughlin ^7j to a multiphase case. The subsequent development of the Lax-Levermore 
theory JH], [HI has shown that the local waveform in the oscillating regions is indeed 
described by the multiphase solutions. The weak limits in the Lax-Levermore problem then 
can be regarded as the result of the Whitham averaging over these solutions. 

A direct formulation of the dispersive shock problem in terms of the Whitham equations 
had been proposed (although without rigorous justification) much earlier by Gurevich and 
Pitaevskii (GP) ^H] who supplied the Whitham system for the KdV equation with natural 
matching conditions at the boundaries separating smooth and oscillating regions and solved 
the problem analytically for the initial data in the form of a step. The discovery by Tsarev 
of the generalised hodograph method [201, an d Krichever's algebro-geometric construction 
[21J allowed for obtaining of many new exact global solutions to the Whitham equations (see 
for instance [22], !2H] and references therein). Availability of such solutions heavily relies on 
the integrability of the Whitham systems considered and, in particular, on the existence of 
the Riemann invariant representation. 

Thus, either way, integrability seems to be an essential part of the analytic theory of 
dispersive shocks. However, in many physically relevant situations, integrable systems al- 
though providing valuable insight into the qualitative properties of nonlinear dispersive wave 
propagation fail to yield satisfactory quantitative agreement with experimental data (see for 
instance discussion in [25 i n the context of finite-amplitude internal shallow water waves). 
This explains the growing interest in the derivation of relatively simple non-integrable mod- 
els which, while providing a more accurate description of physical effects, still seem to be 
amenable to analytic treatment. In the context of the dispersive shock theory this interest 
is supported by a strong numerical evidence that such features of zero-dispersion limits of 
integrable systems as formation of oscillatory zones and weak convergence are true for the 
systems that are not completely integrable but are structurally similar to integrable ones 
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(see for instance ^3] [21], GUI)- On the other hand, the Whitham method used in the direct 
Gurevich-Pitaevskii description of dispersive shocks does not require integrability from the 
governing system. What is needed is just the existence of the periodic travelling wave solu- 
tions and availability of sufficient number of conservation laws. So, in the absence of rigorous 
general approach it seems a natural idea to assume the single-phase Whitham description for 
a "non-integrable" dispersive shock and explore analytic consequences of such an assump- 
tion. Hyperbolicity of the Whitham equations (which also has to be assumed, for instance, 
on the grounds of numerical evidence of modulational stability) would be an essential part 
of such a problem formulation. Of course, such a heuristic approach requires validation. 
This can be made by (i) testing its results on integrable systems where exact solutions are 
available and (ii) by comparison with available numerical results for non-integrable systems. 

In this paper, we develop the outlined approach for the the decay of an initial step 
problem for a broad class of systems that one may characterise as strictly hyperbolic modified 
by weak dispersion. It is clear that in the step decomposition problem the solution of 
the (quasilinear) Whitham equations must depend on x/t alone, which implies principal 
availability of several integrals of motion for their self-similar reductions. It is also clear that 
existence of the similarity solutions does not rely on complete integrability or existence of the 
Riemann invariants for the Whitham equations ( this, unfortunately, does not guarantee that 
such a solution would be available analytically in the absence of the Riemann invariants). 
However, the Whitham equations being obtained by averaging over periodic family have a 
number of special properties distinguishing them from the general class of hyperbolic quasi- 
linear systems. In particular, they allow for exact reductions to the original dispersionless 
equations in the zero-amplitude (linear) and zero-wavelength (soliton) limits. This general 
property together with the self-similarity imposes a number of restrictions on possible values 
of the modulation parameters at the linear and soliton edges of the dispersive shock defined 
by natural boundary conditions. We derive these restrictions by finding a set of integrals 
along the edge characteristics of the Whitham expansion fan. These integrals, of course, are 
equivalent to the basic integrals for the similarity solution evaluated at the edge points. A 
surprising fact is that, the edge parameters are in a simple and universal way expressed in 
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terms of the linear dispersion relation of the system and one of its nonlinear "dispersionless" 
characteristic velocities. 

The obtained set of restrictions can be viewed as a "dispersive" replacement for the 
classical shock conditions and includes for a general bi-directional case 

(i) a condition for admissible jumps for hydrodynamic variables across the dispersive 
shock (which does not coincide with the classical jump condition); 

(ii) the speeds of the dispersive shock edges as functions of a given jump satisfying (i); 
and, since the parameters in (ii) turn out to be determined not uniquely, 

(hi) the set of inequalities selecting the unique valid set of parameters (ii) and providing 
consistency of the whole construction. These inequalities represent an analog of entropy 
conditions known in classical gas dynamics. 

Thus, the obtained conditions allow one to derive main quantitative characteristics of the 
dispersive shock transition bypassing the integration of the Whitham system. In particular, 
they allow for obtaining the amplitude of the leading (or trailing depending on the sign of 
the dispersion in the system) solitary wave in the dispersive shock, the major parameter in 
observational/experimental data. For instance, the conditions readily yield the well-known 
result of the original Gurevich-Pitaevskii paper which reads that the amplitude of the 
lead solitary wave in the KdV dispersive shock is 2A where A is the value of the initial step. 

We apply the obtained general dispersive shock conditions to two integrable (KdV and 
defocusing mKdV) equations and one non-integrable system describing fully nonlinear ion- 
acoustic waves in collisionless plasma. In all cases our description is in complete agreement 
with previous analytic/numerical results. We also discuss the accuracy and some restrictions 
of the developed approach. 

2 Gurevich-Pitaevskii problem for the KdV equation 

We start with an exposition of the theory of Gurevich and Pitaevskii (GP) ^H] formulated 
originally for the KdV equation and later generalised to other integrable equations by differ- 
ent authors (see j2Z| for a detailed introduction and many useful references). Although all 
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formulas in this section are known very well, it is instructive for the purposes of this paper to 
have them handy as they will allow us later to draw parallels with "non-integrable" theory. 
We take the KdV equation in the form 

u t + uu x + u xxx = . (1) 

Let the initial perturbation u(x, 0) = Uq(x) have the form of a large-scale (Ax ^> 1) mono- 
tonically decreasing function with a single inflection point. As a special important case, one 
considers a smooth step 

m (— oo) = u~ , u (+oo) = u + , u' (x) < , (2) 

so that the characteristic width of the transition region, say / is finite. Of course, without 
loss of generality the initial step for the KdV equation can be normalized to unity but for 
our purposes it is convenient to retain arbitrary values for u~,u + ; we only assume u~ > u + . 

The qualitative picture of the KdV evolution of a smooth monotonic profile is as follows. 
During the initial stage of evolution, \u xxx \ <C \uu x \ and one can neglect the dispersive term 
in the KdV equation (0). The evolution at this stage is approximately described by the 
dispersionless (classical) limit of the KdV equation: 

u « (3(x, t): & + f3f3 x = , p(x, 0) = u Q (x) . (3) 

The evolution (J3J) leads to a gradient catastrophe : t — > tf, : j3 x — ► — oo. For t > tb 
(without loss of generality we put t b = 0), the dispersive term u xxx should be taken into 
account in the vicinity of the breaking point and, as a result, the regularization of the 
singularity happens through the generation of small-scale nonlinear oscillations confined to 
a finite, albeit expanding, space region. This oscillatory structure represents a dispersive 
analog of a shock wave. Near the leading edge of such a dispersive shock (or an undular 
bore in a different terminology) the oscillations are close to successive solitary waves and in 
the vicinity of the trailing edge they are nearly linear (see Fig. 1). This structure has been 
recovered in a rigorous theory using the complete integrability of the KdV equation (see 
for instance [21]) but in the Gurevich-Pitaevskii approach it is an assumption based, say, 
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on the results of numerical simulations. This approach is consistent with the aims of this 
paper which develops a way to deal with generally non-integrable systems, where the 1ST 
formalism is not available in principle. 




Figure 1: Oscillatory structure of the dispersive shock evolving from the initial step (dashed 
line) . 

We shall model the local waveform of the dispersive shock by the single-phase periodic 
solution of the KdV equation travelling with constant velocity c: u(x,t) = u(9), 9 = x — ct. 
This solution is specified by the ordinary differential equation 



(u e ) 2 = -G{u) , u{9 + 2n/k) = u{9) , 



(4) 



where 



G(u) = i(w - ui}(u - u 2 ){u - u 3 ) , 
o 



(5) 



U3 > u 2 > ui being constants of integration. The phase velocity c and the wavenumber k 
are expressed in terms of the roots Uj as 

-1 

(6) 
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1 / f du \ 71 (U 3 -Ml) 1/2 

c=3(« 1+ » 2+ »3), ^l/^= =^ K(m) . 

V2 / 

where K(m) is the complete elliptic integral of the first kind. The modulus m and the 
amplitude a of the travelling wave are expressed in terms of {uj} as 



m 



u 3 - u 2 
u 3 - Ui 



a = u 3 — U2- 



(7) 



Eq. (jlj) is integrated in terms of Jacobian elliptic function cn(£; m) to give, up to an arbitrary 
phase shift, 

u(x, t) = u 2 + a cn 2 ^^/ ^ g ^ ( x ~ c ^)i m ^ ■ (8) 
When m — > (tt2 — » U3) the solution (|SJ) turns into the vanishing amplitude harmonic wave 
u(x, t) ~ u 3 — a sm 2 [k (x — cot)] , a = u 3 — w 2 <C 1 , (9) 

where fc = fc(tti, 1*3, 1*3), c = c(ui,u 3 ,u 3 ). The relationship between c and k is obtained 
from Eqs. © considered in the limit U2 — > M3: 

c = m 3 - (10) 

and agrees with the KdV linear dispersion relation for small-amplitude waves propagating 
on the background u = u 3 . 

When m — > 1 (1*2 — > Wi), the cnoidal wave (JHJ) turns into a soliton 

u s (x,t) —Ui + a s sech 2 [A/a s /3(a; - c s t)) , (11) 

which speed c s = c(ux, Ui, u 3 ) is connected with the amplitude a s by 

c s = Mi + a s /3. (12) 

In Eqs. dH]) , (fT^j) . ui plays the role of a background which can be put equal to zero for the 
uniform solutions (as well as u 3 in (JH1), (jl()|) ). However, this is not the case if one is interested 
in modulated travelling waves where the local integrals Uj become slow functions of x and 
t and the background term can not be eliminated by the passage to the moving reference 
frame. This simple fact will play an important role in what follows. 

We now consider slowly modulated cnoidal waves by letting the constants of integration 
Uj be functions of x and t on a large spatio-temporal scale Ax 3> 1, At ^> 1. Then the 
evolution equations for Uj (x, t) can be obtained by averaging any three KdV conservation 
laws d t Pj + d x Qj = over the period of the travelling wave (jHJ) and are referred to as 
the modulation, or Whitham, equations. The averaging is made according to (UJ) as 

2n/k u 3 

k f k [ F(u) 

F(u h u 2 ,u 3 ) = — / F(u(9\u 1 ,U2,u 3 ))d9 = - —===du. (13) 
2vr J TV J J-G(u) 



u 2 



In particular, the mean is calculated as 

u = ui + 2(u 3 - ui)E(m) /K(m) , (14) 

where E(m) is the complete elliptic integral of the second kind. 

As a result, the KdV modulation system is obtained in a conservative form 

d— d — 

— P j (u 1 ,u 2 ,u 3 ) + —Q j (ui,u 2 ,u 3 ) = 0, j = 1,2,3, (15) 

where Pj, Qj are expressed in terms of the complete elliptic integrals. 

One of the modulation equations can be replaced by the so-called wave number conser- 
vation law, which represents a consistency condition in the formal perturbation procedure 
equivalent to the Whitham averaging (see e.g. [2H]), 

a + &= ' (16) 

where u> = kc is the frequency. Of course, Eq. (fTBj) can be obtained as a consequence of the 
three modulation equations (fTo")) fTo] . 

It has been discovered in ^H] that, upon introducing symmetric combinations 

a «l + «2 a U 1 +U 3 U 2 + U 3 

Pi - — 2 — ' ^ 2 ~~ — 2 ' ~~ — 2 ^ ' 

the system ()15|) reduces to its diagonal (Riemann) form 

^+VS-(A,fc,&)|| = l J = 1,2,3. (18) 

where the characteristic velocities V 3 > V 2 > V\ are certain combinations of complete el- 
liptic integrals of the first and the second kind. They can be conveniently represented in a 
"potential" form |U|, [H] 

dM dk_ 

following from the consistency of wave number conservation law (J16|) with the Riemann 
system ()18j) . Here k, m and c are expressed in terms of the Riemann invariants as 

yfQ K(m) P^- Pi 3 
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Of course, the existence of the Riemann invariants for the quasi-linear system of the third 
order (fTK|) is a nontrivial fact and this feature is connected with the preservation of the KdV 
integrability under the averaging. A regular way of obtaining the KdV-Whitham system in 
the Riemann form has been developed by Flaschka, Forest and McLaughlin J7| using the 
methods of finite-gap integration. For the single-phase case of our interest an elementary 
general approach not requiring algebraic geometry can be found in the monograph ^ZT\ . The 
explicit expressions for Vj in terms of the complete elliptic integrals of the first and the 
second kind can be found elsewhere (see for instance |32j, |2H])- We will present here only 
their limiting properties. 
In the harmonic limit 

m = 0: V 3 = (3 3 , V 2 = V l =2f3 l -fo (21) 

so that the Whitham system reduces to 

In the soliton limit 

m = l: V 2 = V 3 = ^(p 1 + 2fc), V l =(3 l (23) 
and the Whitham system reduces to 

Thus, the Whitham system admits nontrivial exact reductions to the dispersionless limit via 
singular limiting transitions fa — > fa (linear limit) and fa —>■ fa (soliton limit). In both 
limits one of the Whitham equations converts into the Hopf equation, while the remaining 
two merge into one for the Riemann invariant along a double characteristics. Such a special 
structure of the Whitham equations makes it possible to formulate the following natural 
matching problem for (J18j) |19j . 

Let the upper (x,t) half-plane be split into three domains : {(x, t > 0) : (— oc,x~(t)) U 
\x~(t), x + (t)} U (x + (t), +oo)} (see Fig. 2), in which the solution is governed by different 
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equations: outside the interval [x~ (t) , x + (t)} it is governed by the dispersionless limit (j3J) 
while within the interval [x~(t), x + (t)] the dynamics is described by the Whitham equations 
(j!8|) so that the following matching conditions must be satisfied: 

x = x-{t): /3 2 = Pi, fo = /3 
x = x + (t): P 2 = fa, Pi = (3 

where ft(x,t) is the solution of the Hopf equation © and the (free) boundaries x ± {t) are 
unknown at the onset. 




Figure 2: Splitting of the (x, t) - plane in the Gurevich - Pitaevskii problem 

Remark. One can notice that formulated in this way, the GP problem contains an 
implicit assumption about the spatial structure of the dispersive shock, namely, it identifies 
the leading edge with the soliton and the trailing edge with the linear wave from the very 
beginning. Of course, for the KdV equation this wave pattern can be inferred from previous 
numerical simulations or from simple physical reasoning that owing to the negative sign of 
the dispersion in the KdV equation, the longer waves must propagate with greater speed. 
However, from mathematical point of view, such an assumption must be confirmed or rejected 
by the actual solution so, strictly speaking, at this point we actually don't know if x + > x~ . 

For the initial data with a single breaking point, the conditions (|25)l uniquely define the 
global solution ftj(x, t) of the Whitham equations (fTHj) (see |2H])- It follows from Eq. (|23j) and 
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the limiting properties of the Whitham velocities (|2ip. ()23|) that the free boundaries x ± (t) 
are determined by the double eigenvalues of the Whitham system for m = (x = x~(t)) 
and m = 1 (x = x + (t)) and are found from the ordinary differential equations 

dx-/dt = V 2 ((3 u (3 u (3 3 )\ x=x -=V l ((3 u (3 u (3 3 )\ x=x - = V-(x-,t), (26) 



dx + /dt = V 2 ({3 u f3 3 ,{3 3 )\ x=x+ = V 3 ({3 U (3 3 ,{3 3 )\ X=X+ = V + (x + ,t) (27) 

where {[3j = j3j(x,t)} is the solution of the GP problem. We note that the curves x ± (t) are 
sometimes referred to as the phase transition boundaries since they separate the regions of 
the zero-phase (external smooth flow) and the single phase (dispersive shock) solutions. 

The local integrability of the Riemann system ()18j) is based on certain relationships be- 
tween the characteristic velocities and has been established by Tsarev [20] who proposed 
a generalisation of the classical hodograph method applicable to the hydro dynamic-type 
systems with the number of field variables exceeding two. There is, however, a special 
important case when the generalised hodograph transform degenerates and the Tsarev in- 
tegrability scheme becomes redundant. This is the case of the decay of an initial step 
which implies that the modulation variables are the functions of s = x/t alone for t » I, 
where / is the characteristic width of the (smooth) initial step. In this case, the GP problem 
has the asymptotic solution in the form of the centred simple wave in which all but one 
Riemann invariants are constant: 

Pi = u + , #, = iT, V 2 {u~,(3 2 ,u + ) = s, (28) 

or explicitly, 

u - + 2u + + mA - — } = - , 29 

3 3 h{m)/ K(m) — (1 — m) t 



where A = u~ — u + is the magnitude of the initial jump. Since the solution ()29|) represents a 
characteristic fan it never breaks for t > and therefore is global. The self-similar boundaries 
s~ and s + of the dispersive shock are found from the solution (|2*9*|) by putting in it m = 
and m = 1 respectively: 

s - = s (0) = u + - A , s + = s(l) = u + + ^A , (30) 
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One can see that s + > s so the assumed wave pattern with the leading soliton and trailing 
linear wavepacket was indeed correct. Now the amplitude of the lead soliton is simply 

a+ = 2(/5 3 -/3i)=2A, (31) 

whereas the value of the wavenumber of the trailing wave packet follows from Eq. ([2 Oh) 
evaluated for f3 2 = (3\ which yields 

k ~= \l\w*-h)= vF- (32) 

We see that the simple form (J28|) of the analytic solution to the GP problem is possible 
owing to the availability of the Riemann invariant form (|18j) which admits (3j = constant 
as an exact solution. The formulas (JBTlj) - ()32|) are obtained as consequences of this global 
solution. So, since the existence of the Riemann invariants is due to complete integrability 
of the original KdV equation one may conclude that the possibility of obtaining the "global" 
formulas (j3*Uj) - (j3*2^) relies on the integrability as well. This is not in fact so. In the next 
section, we will show that these formulas can be drawn directly from the Whitham system 
in "physical" variables using some very general properties of the Whitham systems subjected 
to natural boundary (matching) conditions of the GP type. Actually, Eqs. (p?Tj) . ()32|) will 
be obtained without making use of the Riemann invariant form for the Whitham equations, 
and bypassing their explicit integration. 

3 "Non-integrable" reformulation of the Gurevich - 
Pitaevskii problem 

3.1 Formulation 

Let us reformulate the GP problem in terms of the Whitham system in its original, non- 
diagonal form (|15|) . Along with the set of natural "mathematical" variables u±, u 2 , u% we will 
use an equivalent set of "physical" variables u,k,a which are connected with ui,u 2 ,U3 by 
means of Eqs. 0, ©, ©. 
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From here on we are not going to exploit subtle algebraic properties underlying the inte- 
grability of Instead, we will use some general properties of the system (jTHJ) connected 
with its "averaged" origin and distinguishing it from a general class of hyperbolic quasi-linear 
systems of the third order. In all relevant cases we will indicate if the obtained relationships 
could be extended to a more general context. 

We start with the KdV modulation system in conservative form (see JT^J for instance) 



du d I u 2 \ „ d ( u 2 \ d ( u 3 u 2 n \ „ dk du>(u, k, a) 

(33) 

It is not difficult to show that the averaging (J13|) over the period of the travelling wave 
specified by Eqs. (HJ), (JSJ) implies the following general relationships for the averaged variables 
in the harmonic (m = 0) and soliton (m = 1) limits 



F(u)\ U2=U3 = F(u 3 ); F(u)\ U2=ul = F( Ul ); (34) 

In particular, this implies that u{ui, u 3 , w 3 ) = u 3 and u(ui, Ui, u 3 ) = u\ (this can also be 
directly seen from the explicit expression (j!4j) ). Then, since = u 3 is equivalent to a = 
and U2 = u\ to k = 0, the relationships (pMj) assume in the physical variables the form 



F(u)\ a=0 = F(u); F(u)\ k=0 = F(u); (35) 

The relationships ()35|) immediately imply that in both harmonic and soliton limits, u 2 = u 2 
hence the first modulation equation (J33j) turns into the Hopf equation for the mean value 
u t + uu x = 0. It can also be readily show shown using (jlj), (JT3J) that the same is true for the 
second modulation equation (|33J) as well. Therefore: 

i) the modulation system (J33|) admits exact reduction to a lower order system both for 
a = and k = 0. The limiting transitions a — >• and k — > are, therefore, singular; 

ii) the limiting transitions a — * and k — > provide two possible ways of passage to the 
dispersionless limit in the modulation equations. 

We note that the same conclusions have been inferred from the Riemann equations (|22j) . 
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Of course, the described limiting behaviour is quite expected for a modulation system 
since in both afore-mentioned limits the oscillations do not contribute into the averaging 
hence the averaged (semi-classical) system must agree with non-oscillatory classical limit 
of the original equation. It is also clear that this fact is not unique to the KdV equation 
and must hold for other dispersive-hydrodynamics systems supporting the periodic travelling 
waves with the quadratic behaviour of the potential curve G(u) © near its extrema. One 
can see, in particular, that this (generic for weakly dispersive systems) behaviour of the 
potential curve guarantees exponential decay for the soliton solutions so that they do not 
contribute into the averaging. 

Now the GP natural boundary conditions (J25)l can be reformulated in terms of the match- 
ing of the mean flow in the oscillatory region with the smooth external flow 

x = x~(t) : a = 0, u = 8 , 

(36) 

x = x + (t) : k — 0, u — 8, 

where (3(x, t) is defined by the classical limit (j3J). Using explicit expression for u ((Hj) and the 
relationships (fTTj) one can easily see that the conditions (f3T)|) indeed equivalent to (J25J) . For 
the particular case of the decay of an initial discontinuity ((21), the conditions assume 
the form 

x = s~t : a = , u = u~ , 

(37) 

x = s + t : k — , u = u + . 
We note that when passing from ((36|) to (|37j) we have made an implicit assumption about 
the structure of the dispersive shock by implying that s + > s~ . 

The boundaries x ± {t) of the modulation solution are defined by Eqs. (J27J) in terms 
of the characteristic velocities of the modulation systems where either m = (trailing edge) 
or m = 1 (leading edge). This definition, of course, does not rely on the existence of the 
Riemann invariants and must hold in general case. This can be explained in the following way. 
It is clear that in order to provide continuous matching of the solutions of two (consistent) 
quasilinear hyperbolic systems of different order (the Whitham system and the Hopf equation 
in our case) the matching lines must necessarily be the multiple characteristics for the system 
with a higher order. Thus, the natural matching conditions (|3l)j) must be supplemented with 

16 



the definition of the boundaries x ± (t) in terms of the double eigenvalues of the modulation 
system (see [3B] for a detailed description of the characteristics behaviour in the GP problem). 
We represent the Whitham system (JT5j) in a generic quasi-linear form 

y t + A(y)y x = 0, (38) 

where y = (u, k, a) T and A(y) is the coefficient matrix. As is well known (see for instance 
|34j . [35] ) the quasilinear system (}38j) can be represented in a characteristic form 

V/,W^-0 ^ = A + t/A m -i23 (39) 
i=i 

where the characteristic velocities V}(y) are the eigenvalues of the matrix A(y), and b^ m ^(y) 
is its left eigenvector corresponding to m-th eigenvalue: h^A = V m h( m \ We assume that 
the characteristic velocities Vi(y) < V 2 (y) < V^y) are real for the solutions of our interest 
so the system ()38j) is hyperbolic (of course, for the KdV equation, the hyperbolicity of 
the Whitham system is a proven fact |3E] but in our, "non-integrable" approach it is an 
assumption) . 

It can be readily seen using an equivalent characteristic representation ()39|) in "mathe- 
matical" coordinates, y i— > y* = (^1,^2,^3), that the fact that the modulation system (fHKJ) 
admits exact reductions for 112 = u% (a = 0) and 112 — U\ (k — 0)implies that the multiple 
roots of the potential curve G{u) necessarily correspond to multiple eigenvalues of the av- 
eraged system. This, again, is a consequence of the properties of the Whitham averaging 
(j!3)) . It follows from the characteristic velocities ordering that the double eigenvalues can 
be either V2 = V3 or Vi = V\. The remaining single eigenvalue in both cases is u since the 
Hopf equation is an exact reduction of the Whitham system in both limits. If the explicit 
expressions for the characteristic velocities are known, the correspondence between the dou- 
ble roots of the travelling wave potential function G(u) and the double eigenvalues of the 
averaged system is established directly (see (j21J) . (J2SI))- Alternatively, one can perform an 
asymptotic analysis of the modulation system for m C 1 (a C 1) and (1 - m) < 1 (fe < 1) 
and establish the sought correspondence. Either way, even without additional analysis it is 
clear that, due to the accepted ordering of Vj's, the double eigenvalues should be identified 
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with the values of the middle characteristic velocity V 2 in the harmonic and the soliton lim- 
its. So we define the boundaries x ± (t) of the KdV modulation solution by the characteristic 
equations (cf. (Poll. (1771)) 

dx~/dt = Viiu, k, 0)\ x=x - , dx + /dt = V^w, 0, a)\ x=x + . (40) 

For the decay of an initial discontinuity Eqs. (f4*U|) . in view of assume the form 

s~ = V 2 (u-,k-,0), s + = V 2 (u + ,0,a + ) (41) 

where the values a + and k~ are to be found from the solution of the GP problem. 

In the subsequent sections we will find simple effective expressions for double characteris- 
tics of the Whitham system using physical variables u,k,a and bypassing the full eigenvalue 
analysis. Moreover, we will show that the edge parameters a + , k~ and hence s + , s~ can be 
found from a certain set of conditions not involving the global integration of the Whitham 
equations. As a result, we will present a way to "fit" the dispersive shock into the solution 
of the dispersionless limit equations without the detailed analysis of its internal oscillatory 
structure in the same manner as the traditional shock is embedded in the solution of the 
inviscid equations of ideal gas. As will follow from the construction, the availability of these 
conditions does not depend on the existence of the Riemann invariants for the Whitham 
system. 

3.2 Trailing edge 

We start with the determination of the trailing edge speed s~ for the problem of the decay of 
an initial discontinuity (J2J) in the KdV equation using the reformulation of the GP problem 
in physical variables made in the previous subsection. The trailing edge is defined for the 
solution of the GP problem by the characteristic where a = 0. Our analysis below will be 
based on the special features of the Cauchy data prescription on characteristics (see e.g. 

It follows from the characteristic form that the differentials of the modulation vari- 
ables along the z-th characteristic are not independent but related by means of the expression 
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b^i (u, k, a)diU + b% {u, k, a)dik + fog (u, k, a)dia = . (42) 

Let a = on this characteristic (which necessarily implies that it is a double characteristic 
- see explanation in the previous section). Since a = must be an exact solution of the 
modulation system, the ordinary differential equation 

bf (u, k, 0)diU + bf {u, k, 0)d { k = (43) 

necessarily represents a characteristic equation for the reduction as a = of the full mod- 
ulation system. This equation can always be integrated using an integrating factor to give 
the relationship between admissible values of u and k on the characteristic, i.e: 

Mu,k)=C , (44) 

Co being constant of integration. The point is that this relationship is local, so it does not 
depend on the specific solution under study but is determined (up to a constant) only by 
the coefficients of the Whitham system evaluated for a = 0. Therefore, the integral (J44|) 
can be derived, bypassing the technically involved route via full characteristic form ()42|). by 
a direct substitution of the ansatz k = k(u) into the reduction of the Whitham system for 
a = 0. We also remark that the value Cq represents, in fact, a Riemann invariant for this 
2x2 reduction. Of course, the existence of this Riemann invariant does not depend on the 
existence of the Riemann invariants for the full Whitham system fTo*j) . 

It follows from the relationship (jSSK) and Eqs. (jHJ), (fT4^) considered for u 2 = u 3 , that 

lim u 2 = u 2 , lima; = lim ck = ku — k 3 = uj${u, k) . (45) 

a— >0 a— >0 U2— >«3 

Obviously, the second expression in Eq. (j43|) represents the KdV linear dispersion relation, 
where the linearization is made about the slowly varying mean background u(x,t). Now we 
can immediately obtain the reduction of the Whitham equations (|3*3*j) 

a = 0, u t + uu x = 0, k t + (u (u, k)) x = . (46) 

The equations (J46|) form a closed system. We emphasize that the system (J46j) is not an 
asymptotic system, it is an exact reduction of the full modulation system. Of course, the 
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reduction ()46|) is equivalent to the limiting diagonal system (}2*2*|) (which is readily verified by 
a direct calculation) but here it was derived directly from the original modulation equations 
()33|) without using the Riemann invariant form. Actually, the limit ()46|) has a clear meaning 
of the formal modulation system for a zero-amplitude wave packet propagating on a slowly 
varying background flow u and can be postulated on the physical level of reasoning directly, 
even without turning to the Whitham equations (J3*3*j) and the formal derivation of the linear 
dispersion relation ()45|) for modulated waves via nonlinear travelling wave solution. Now, 
looking for the integral k(u) of Eqs. (f4T)J) we arrive at the ordinary differential equation 
(which is an equivalent of Eq. (|IHJ)), 

dk Ouiq/Ou 
du u — duj 1 dk ' 

which, upon substituting the linear dispersion relation from (|43|) is readily integrated to give 

k=^(u + C 1 ). (48) 

Thus Eq. ()48)1 is a relationship between values of k and u on the characteristic of the Whitham 
system on which a = (there is no relation to any particular global solution yet). 

Now we apply the relationship (}4*5|) to the GP problem (}3Tj) . where the characteristic 
where a = is associated with the trailing edge x = s~t. First, we find the constant C\ from 
the second boundary condition (|37|) which, in the space of the field variables prescribes k = 
when u = u + (this condition does not contain the amplitude a so it must hold everywhere in 
the plane k = of the three-dimensional space with the coordinates {u, k, a}, including the 
line (u, 0, 0) as well). Thus C\ = —u + . Then, putting u = u~ in (j4~%|) we obtain the value of 
the wavenumber at the trailing edge of the dispersive shock x = s~t, 

k~ = V2A/3, (49) 

where A = u~ — u + is the jump across the dispersive shock. This agrees with the value of 
k~ given by Eq. (J32j) obtained as a consequence of the full modulation solution (J61)) . 

Now we find the self-similar co-ordinate (the speed) of the trailing edge s~ which, ac- 
cording to (j4*Uj) . is calculated as the double characteristic velocity of the Whitham system 
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for the solution of the GP problem. It readily follows from the limiting system (|46jh that its 
characteristic velocities are 

u, ^r(u,k) = u-3k 2 . (50) 
ok 

The latter, of course, is the group velocity of the linear wave packet propagating on the 
varying mean flow background u(x,t), which is perfectly reasonable from physical point of 
view. Now we need to identify one of the velocities (joTlj) with the double eigenvalue V^m, k, 0) 
in ([4 lb) to evaluate the trailing edge speed. It is clear from physical reasoning that this must 
be the linear group velocity. Later we will present additional formal conditions eliminating 
any possible ambiguity in the identification of the double eigenvalue. Of course, for the 
KdV case it is not difficult to show using explicit formulas (JHJ), © that expressions (jBTil) 
are equivalent to those obtained from the Riemann form of the Whitham equations in the 
harmonic limit (see Eq. (J2T|) ) and our identification V^w, k, 0) = duo/dk is indeed correct. 
Thus, we get from Eq. (|41b) for the trailing edge 

s " = w ( ""*"»- (51) 

In this form the equation of the trailing edge can be interpreted as a kinematic bound- 
ary condition for linear modulated waves and can be postulated as a part of the problem 
formulation. 

Setting (JSIJ), i 11 ^ (JHH) we obtain for the trailing edge 

s~ = u~ - 3(&T) 2 = u + - A. (52) 
This value, again, agrees with the full solution of the GP problem (see Eq. ([3Uj)). 

3.3 Leading edge 

In principle, one could proceed with the leading edge in the same fashion as we did with the 
trailing edge i.e. by indirect integration of the full characteristic 1-form along a characteristic 
where k = 0, where this form degenerates into 

(u, 0, a)diU + 62 (^j 0> a)did = . (53) 
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That would imply finding the integral a(u) for the exact zero-wavenumber reduction of the 
modulation system. Such a reduction, however, must include full modulation equation for the 
soliton amplitude, which can be obtained without too much trouble for the KdV equation 
but is not readily available for more complicated systems, especially if the solitary wave 
solution can not be found explicitly One should remark that the "obvious" universal soliton 
amplitude equation a t + c s (a)a x = (see [31], Sec. 16.6) is correct only if the background 
flow u is constant. The full amplitude equation for solitary waves (see [37], [HE] for instance) 
takes into account variations in all modulation parameters and inevitably contains the term 
proportional to u x which is vital for deriving the characteristic equation (j53|) . As a result, 
unlike the wave number conservation law, the full soliton amplitude equation does not have an 
explicit universal form and a straightforward application of the method proposed in Section 
3.2 would be essentially equivalent to the consideration of the general characteristic equation 
(j42j) in the soliton limit. For "real" non-integrable systems that could be a very difficult 
technical task. This complication, however, can be bypassed by introducing a conjugate 
wavenumber defined for the KdV as 
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t( Ul ,« a , Ua )=H l-^j R(ml) , m=l-m (54) 

- instead of the amplitude a and the ratio A = k/k instead of the original wavenumber k 
©• The new set of modulation variables we are going to use is (u, A,k). It is readily seen 
that the soliton limit k — (to = 1) corresponds to A = 0, k = a s /3 = 0(1) in new 
variables. On the other hand, in the harmonic limit a = (m = 0) one has k = 0. So k 
indeed plays the role analogous to the amplitude. For further convenience, we rewrite the 
matching conditions (|H7jl using new variables 

x = s~t : k — , u = u~ , 

(55) 

x = s + t : A = , u = u + . 

Now, as was shown in Section 3.1, the two first averaged conservation laws (j33|l in the soliton 
limit reduce to the Hopf equation: 

A = : u t + uu x = 0. (56) 
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To obtain the equation for k in the soliton limit we set the new variables in the wave number 
conservation law in to obtain: 

kA t + uA x + A(k t + uj x ) = 0, (57) 

where Co = u(u, A, k) = kc is the conjugate frequency. 

Now, using the arguments identical to those in the previous subsection we infer that 
if A = on some characteristic then the values of k and u on this characteristic must be 
connected by a local relationship k(u) which must not depend on the particular solution. 
We consider equation (JSTj) in the small vicinity of the mentioned characteristic where A < 1. 
Then, using k = k(u) and Eq. (}56| for the leading order, we obtain an asymptotic equation 

OA u s dA Adu\dkfdu s \ du s \ /n / A dA\ , , 



dt k dx k dx \ du \ dk J du \ \ dx J 

where the reduction uj s (u,k) = u)(u,0,k) can be called a soliton dispersion relation. Since 
the sought characteristic integral k{u) must not depend on the way A tends to zero, the 
expression in the brackets (which does not depend on A) must be identically zero and Eq. (j55jl 
splits into 

dk du s /du , . 

— = = 59 

du u- dus/dk 

and 

8A + a.aA = £ ,/ A m_ 



dt k dx \ dx J 
Eqs. (|59|l. (jnOj) have been obtained in [39J using somewhat more particular arguments. One 
can not help noticing that the equations (|37j) and ()59|) for the characteristic integrals k(u) 
(linear characteristic) and k(u) (soliton characteristic) are identical in terms of the corre- 
sponding dispersion relations u (u, k) and Co s (u, k). The latter, however, is yet to be found. 

It follows from (jfiOjl . (|5()jl that in the soliton limit the characteristic velocities of the 
Whitham system are 

J, «• (61) 

One can see from the definition of Ca s that the characteristic velocity u s /k coincides with 
the soliton speed c s (u, k) = c(ux,Ui,u 3 ), which is, of course, expected. We identify c s with 
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the double characteristic velocity V 2 = V 3 of the full Whitham system in the soliton limit, 
then the classical speed u would correspond to V\. The conditions eliminating possible 
ambiguity in the characteristic velocity identification for the reduced Whitham system will 
be presented later. Now, similarly to the trailing edge case, the definition of the leading edge 
in Eq. (J41j) can be represented as a kinematic boundary condition for the lead soliton (cf. 

m ) 

s+ = c s (u + ,k + ), (62) 

where k + = k(u + ). One can see that in this form the definition of the leading edge can be 
adopted part of the problem formulation. 

At this point, one may take advantage of the explicit expressions (JHj), © for u(u\, u 2 , u 3 ), 
k(ui, u 2 , u 3 ) and c(ui, u 2 , u 3 ) to find u a (u, k) and solve Eq. (|5^j) using the matching conditions 
(|55|l . Such an attractive route, however, might not be readily (if at all) available for actual 
non-integrable systems, where the physical modulation parameters often can not be expressed 
in a simple way in terms of the roots of the potential curve. So we proceed with some 
apparently "roundabout" way, which would be more universally applicable to other systems. 

We observe that k and Cj = kc defined by (|H4"j) and (JBk) as functions of the roots ui,u 2 , u 3 
can be viewed as the wavenumber and the frequency in the conjugate travelling wave associ- 
ated with the same set of roots Uj as in Eq. but inverted potential curve (now oscillations 
occur between the roots u\ and u 2 ): 

{u § ) 2 = G{u) , u{9 + 2Tx/k) = u{6) . (63) 

where u is new field variable and 6 = x — ct is a new travelling phase associated with "old" 
phase velocity c = {u\ + u 2 + u 3 )/?>. Eq. (|6l3|) can be obtained from (jlj) by the change of 
variables u 1— > u, x 1— > ix, 1 1— > it, which corresponds to a mere change of the dispersion sign 
in the KdV equation (JTJ), 

u- t + uu x - Uxxx = . (64) 

As a matter of fact, in our particular KdV case the functions u{9) and iu(—i8) defined for the 
same set of roots Uj represent the same analytic (elliptic) function in the complex x-plane 
with the periods 27i/k and 2iri/k along the real and the imaginary axes. Generally, of course, 
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this is not the case, so we are not going to take advantage of the analytic properties of u(8) 
here. Instead, we observe that the soliton limit for the original travelling wave equation (j3J 
(112 — > u\) corresponds to the harmonic limit for the conjugate equation ()63|) . Therefore, the 
relation between u s and k can be obtained as a linear dispersion relation for the conjugate 
KdV equation (JMJ) . i.e. 

u s = k < u > +k 3 , (65) 
where the brackets < ... > denote the averaging over the conjugate family (JfiHj) (cf. (fT3]l) 



~ U2 



k f F(u) 

<F> {u h u 2 , us) = - / -y=Ldu . (66) 
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It is not difficult to show by a direct calculation that < F > (u\,ui, U3) = F{u\) . Then, it 
follows from Eq. (J34)) that for w 2 = U\ (k = 0) we have 



<F(u) > \ k=0 = F(u)\ k=0 (67) 

and, in particular, < u >= u. Again, similarly to Eqs. (jHU), (|3*3jl this property is due to 
the quadratic behaviour of the potential curve G{u) ©in the vicinity of the double root 
corresponding to the soliton limit and is not confined to the KdV example alone. 

It follows from Eq. (jBTj) that u s (< u >,k) = u s (u, k) and therefore, the soliton dispersion 
relation can be obtained from the original linear dispersion relation Uq(u, k) by the formal 
change 

k 1— > ik , uo 1— > iuo (68) 
In other words, the soliton dispersion relation is found as 

u s (u, k) = —iuo(u, ik) = ku + k 3 . (69) 

Of course, in view of k(ui, U\, 113) = y/ a s /3 and u{u\,u\,Us) = U\, Eq. (j69|) is equivalent 
to the classic relationship (jllj) between the speed and the amplitude of the KdV soliton 
but here it was obtained in terms of the linear dispersion relation so the outlined procedure 
can be easily generalised to more complicated systems, where the explicit analysis of the 
travelling wave solution is not as readily available as in the KdV case. 
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The possibility of expressing the soliton speed in terms of the linear dispersion relation 
might look somewhat surprising but can be explained by the following simple argument (see 
|40j): the value of u in the soliton tail, which moves with the same speed as the rest of the 
soliton, is very small, so one can, in principle, infer this speed from the linear theory. 

Remark. We emphasize that the obtained relationships between the original and the 
conjugate averaged variables essentially represent algebraic identities between integrals of 
the form (JT5)l and (|66p associated with a given Riemann surface fi 2 = G(u) and do not imply 
any connection between the spatio-temporal modulation dynamics for the original and the 
conjugate equations ((TJ) and (jHlj). 

Now setting Eq. (JBTJj) into Eq. (j59|) we obtain after elementary integration 

k=^(C 2 -u), (70) 

where C\ is constant of integration. We emphasize that Eq. (|7U|) (as well as Eq. (J48|) in the 
linear case ) is a general relationship between values of k and u on the characteristic of the 
Whitham system where A = and is not tied to any particular global solution. 

Now we apply the relationship (J7U|) to the GP problem (|55|). where the characteristic 
on which A = is associated with the leading edge x = s + t. First, we find the constant 
C*2 from the first boundary condition (jo3j) which prescribes k = when u = u~ in the 
three-dimensional space with the coordinates {u,A,k}. This condition does not contain A 
so it must hold on the line {u, 0,0} as well and, therefore, can be applied to (J70|) . Thus 
C*2 = u~ . Then, putting u = u + in (fTUj) we obtain the value of the conjugate wavenumber 
at the trailing edge of the dispersive shock s = s + : 

k + = y/2A/3. (71) 

We observe that, for the KdV equation k + = k~ (see (J49j) ) . Of course, this is not a general 

relationship although one can expect significant symmetry between the expressions for k~ 
and k + in other dispersion-hydrodynamic systems. 

The speed of the leading edge (jB2I) is now calculated with the aid of the conjugate 
dispersion relation (pj9"j) as 

, uj s (u + , k + ) , 2 . ,„ . 

k+ 3 1 ; 
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which agrees with the global solution of the full Whitham system (see Eq. f|3Uf> ). The 
KdV soliton amplitude is connected with its velocity by the relation p2j) . which in our case 
assumes the form s + = u + +a + /3. Hence the lead soliton amplitude is a + = 3(s + — u + ) = 2 A, 
which, again, agrees with the full modulation solution. 

3.4 "Entropy" conditions. 

One can see that the valid solution to the GP problem with u~ > u + must satisfy the 
inequalities 

s~ < u~ , s + > u + , s + > s~ , (73) 

These inequalities ensure that the "external", Hopf characteristics x = u~t and x = u + t 
starting from the x-axis on either side of the dispersive shock region intersect its edges when 
continued in the direction of increasing t thus transferring the initial data into the dispersive 
shock zone s~t < x < s + t (see Fig. 3). The inequalities f73]l represent a dispersive- 
hydrodynamic analog of the classical gas- dynamic entropy conditions jT]. As a matter of fact, 
these inequalities are redundant when the explicit solution of the GP problem is available 
(one can see that conditions ()73j) are indeed satisfied by the boundaries of the solution ()28j)). 
However, in the absence of the full modulation solution, and without physical assumptions 
about the spatial structure of the dispersive shock there is an ambiguity in the determination 
of the edges using the construction proposed in Sections 3.2, 3.3. This ambiguity has been 
pointed out when we made an identification of the leading and trailing edges with the soliton 
and the harmonic wave in the natural boundary conditions (|3T|) and when we identified linear 
group velocity and the soliton velocity with the double eigenvalues of the full Whitham 
system in respective limits. One may speculate about the possible analog of the entropy for 
dispersive hydrodynamics. An attractive candidate is given by the integral 

X 

Soc J kdx'>0. (74) 

— oo 

For the KdV dispersive shock the positive function k(x, t) is supported on the interval s~t < 
x < s + t and one can see that the integral (J85|) increases when crossing the dispersive shock. 
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Figure 3: Classical (Hopf) characteristics transfer initial data into the Whitham zone 
3.5 Geometric interpretation 

We now place the obtained relationships in a more general context. The quasilinear hyper- 
bolic system of the third order (|38jl admits the centred expansion fan solution 

Fi(u, k, a) = h , F 2 (u, k, a) = I 2 , V k (u, k,a) = x/t , (75) 

Ji, I2 being constants and is one of the characteristic velocities so that the solution 
satisfied the GP matching conditions (|37|) (we know that for the KdV equation this is V2 - 
see Eq. fJ28jl ). The two first expressions in (|75jl define, for given two integral surfaces 

in the space of the field variables u,k,a with k > 0, a > (physical restrictions). Their 
intersection yields the solution curve {T(s; Ii, I2) : u = u(s) ,k = k(s) ,a = a(s)}, where 
the parameter s = x/t changes on the interval [s~ , s + ] so that a(s~) = and k(s + ) = 0. 

We now briefly outline how one proceeds with this geometric construction. First we 
observe that, the endpoints of the solution curve T(s) lie on the lines of intersection of two 
integral surfaces parametrised by Ji and I2 with the coordinate planes a = and k = 0. It 
follows from (|75|l that there are four these lines: Fi }2 (u, k,0) = I\ t2 and F lj2 (w, 0, a) = I\ t2 . 
Since the Whitham system admits exact 2x2 hyperbolic reductions for a = and k — 0, 
the equations of the intersection lines can be obtained directly, as the characteristic integrals 
(Riemann invariants) of the reduced systems, i.e as solutions in the form k(u) and a(u) 
respectively. The constants Ji^ are then expressed in terms of the initial parameters u~, u + 
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by applying the boundary conditions (|37jl to the characteristic integrals as it was done in 
Sections 3.2, 3.3. Then, using the kinematic conditions (jSH), (ffi2~j). one evaluates the speeds 
of the trailing and the leading edges and only after that, the valid pair of the endpoints is 
selected by the "entropy" inequalities (J73J). 

Essentially, this general construction has been realised in Sections 3.2, 3.3 by using the 
arguments from the characteristics theory and introducing a different (conjugate) basis of 
the field variables when considering the leading edge. 

4 Dispersive shocks in simple- wave led dispersive equa- 
tions 

It is clear that the presented construction is not restricted by the KdV equation alone, and 
can be naturally generalised to the "KdV-like" nonlinear dispersive systems possessing the 
basic properties necessary for the Whitham averaging. We represent the governing equation 
in the form 

u t + V(v)u x + K 3 [u] = 0, (76) 

where V(u) is a real function and K% is a real differential operator of the third order with 
respect to spatial or mixed derivatives such that the equation (|75|) has real linear dispersion 
relation to = u (k). We assume the following general properties for the equation (1751) : 

(i) it admits the hyperbolic classical (dispersionless) limit obtained formally by introduc- 
ing the stretched independent variables x' = ex, t' = et and tending e — > while assuming 
finiteness of the derivatives with respect to x' and t', 

u t > + V(u)u x > = . (77) 

In terms of the linear dispersion relation this property implies lvq ~ k for k <C 1 and is 
associated with weakly dispersive waves. 

(ii) it possesses at least two conservation laws; 

(iii) it supports periodic travelling waves parametrised by three independent integrals of 
motion such that the travelling wave solution allows for a harmonic (zero-amplitude) and a 
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solitary wave (zero wavenumber) limits. We will assume the following KdV-like behaviour for 
the "potential" function G(u) in the ordinary differential equation (ite) 2 = —G(u) specifying 
the travelling wave solution (cf. Eq. (J5]l): a) The function G(u) has at least three real zeros 
U3 > U2 > U\ such that the oscillations occur between u 2 and u% (the latter is assumed 
just for sake of defmiteness) ; b) In both nearly linear (u 3 — u 2 )/(us — %) <C 1 and nearly 
soliton (it 2 — Ui)/(u 3 — u\) <C 1 configurations the generic (quadratic) asymptotic behaviour 
is assumed: 

G(u) = (u- u x )(u - u 2 )(u 3 - u)G*(u) , (78) 

where 



G*{u) = 0(1) for 



u — u 2 



< 1 (79) 



U 3 - Mi 

so that the limiting transitions (|34j). ()67|) for the mean values can be easily shown to take 
place; 

(iv) the corresponding Whitham system of the third order (two averaged "hydrodynamic" 
conservation equations plus the wave number conservation law) is hyperbolic for the solutions 
under study. 

The equations of the form (J75j) possessing the mentioned properties (i) - (iv) may be 
characterised as simple-wave led weakly dispersive. We note that all assumptions except (iv) 
can usually be explicitly verified for a specific system. We partially address the hyperbolicity 
issue (iv) in the end of this section. 

We consider initial data in the form of an arbitrary step 

Uq(x) = u~ for x < ; u + for x > , . (80) 

and assume that the dispersive shock transition for t > can be modelled by the centred 
expansion fan solution to the Whitham equations. Now one can see that all arguments 
used in the derivation of the dispersive shock conditions for the KdV equation in a "non- 
integrable" reformulation hold for the case being considered: one just needs to replace the 
Hopf "dispersionless" speed with V(u) and the KdV linear dispersion relation with the 
dispersion relation corresponding to Ij7fij) . So we present only the final relationships without 
making any assumptions about mutual position of the soliton and harmonic edges. 
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First we introduce the linear dispersion relation for Eq. (|76|) by considering an infinites- 
imal perturbation of a mean level u 

u^u + Ul e i{kx ~ ut) , ui < 1 , (81) 

which yields u = coq(u, k) to leading order. Let the dispersive shock transition be confined 
to an interval s~t < x < s + t. We introduce two sets of parameters {k~ , s~; k + , s + }i and 
{k~, s~; k + , s + }2 ■ The set {k~ , s~; k + , s + }i is associated with the negative dispersion wave 
pattern when the soliton appears at the leading edge of the dispersive shock. We define 
this set in the following way. First we find two functions k(u) and k(u) from the ordinary 
differential equations 



dk dujQ I du 



du V(u) — duj /dk 
dk dui s /du 



k(u + ) = 0, (82) 



k(u~)=0. (83) 



du V{u) - du s /dk 

Then the values of the wavenumber at the trailing edge k~ and the lead soliton conjugate 
wavenumber k + are found as k~ = k(u~), k + = k(u + ). By definition, they must be real. 
The speeds of the dispersive shock edges are found from the expressions 

.-=^ ( „- fc -), a + = ^n, ( 84) 

where £u s (u, k) = —iu>a(u,ik). The second set of parameters {k~ , s~; k + , s + }2 associated 
with the reversed, positive dispersion pattern, is obtained from the same system (|52*j) - 
(|84l) in which one replaces "— " with "+" which is equivalent to a respective replacing of 
the parameter subscripts directly in the set {-}i- The valid set is selected by the reality 
condition for the wavenumber and by the "entropy" conditions, 

s~ < V(u~) , s + > V(u + ) , s + > s~ . (85) 

In principle, the solution can switch between the sets {-}i and {-}2 depending on the initial 
conditions. "Switching" of the edge parameter sets implies the reversion of the spatial 
structure of the dispersive shock. Another possibility is that for some domain of the initial 
data (w + , u~) both sets {-}i and {-}2 fail to satisfy either reality or "entropy" condition. This 
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would imply that a genuine global solution to the dispersionless limit equations is available 
for such initial data and no breaking occurs. This solution, of course, is a classical rarefaction 
wave. 

We also briefly address the issue of the modulational stability. Our main mathematical 
assumption in this paper is that of the hyperbolicity of the Whitham system for the solutions 
under study. This ensures global modulational stability of the dispersive shock. For non- 
integrable dynamics, establishing the region of hyperbolicity for the modulation system is 
a separate, often technically involved, problem. However, in a more restricted context of 
the dispersive shock description, some effective necessary conditions of global modulational 
stability can be formulated in the following way. 

We evaluate the the frequency and conjugate frequency at the respective boundaries 
of the dispersive shock transition using the solutions of the ordinary differential equations 
(JE2J), (|S3]). We denote these frequencies as u [u~,u + } and uj s [u~,u + ]. Then, given the 
strict hyperbolicity of the dispersionless limit, the necessary criterion for global modulational 
stability of the dispersive shock is given by the conditions 



which are equivalent to the natural requirement for the speeds s~, s + to be real. The 
conditions (J85j) define the domain D in the initial data plane (u + ,u~) which corresponds to 
modulationally stable solutions. It is not clear if the formulated criterion is the sufficient 
condition so the obtained region D[u + ; u~] could be corrected by some additional conditions. 

We note that, for non-integrable equations some additional restrictions on admissible 
values of initial data can occur due to existence conditions for the single-phase travelling wave 
solutions, for instance, due to the presence of the upper bound for the soliton amplitude. 

Example: Decay of a step problem for the defocusing mKdV equation 

As the next simplest example of an effective construction of the dispersive shock transition 

in the simple-wave led dispersive equation we consider the step resolution problem for the 

defocusing mKdV (mKdV(d)) equation 



Im ojq [u , u + ] = 



Im u s [u , u + ] = , 




u t — u 2 u x + u, 



'XXX 
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u(x, 0) = m_ for x < ; u + for x > . 



(88) 



Equation (J57)l is an exactly integrable equation and can be treated by the 1ST method. It 
belongs to the defocusing NLS hierarchy (self-adjoint spectral operator) and its Whitham 
system is known to be hyperbolic jHj. Actually, the modulation system for mKdV(d) 
equation in Riemann invariants is identical to that for the KdV equation (this fact had been 
established in [31] by a direct calculation before the methods of finite-gap modulation theory 
became available). That, however, does not imply that the physical modulation solutions for 
KdV and mKdV(d) in terms of initial data u + and u~ will be the same or even have same 
properties. We will demonstrate now how major qualitative and quantitative characteristics 
of the modulation solution can be derived from the dispersive shock conditions - (JHU 
without derivation of the Whitham system and analysis of its full solution. 

We infer from flS7|) that V(u) = -u 2 , u Q (u, k) = -ku 2 - k 3 . Then using (JBD - USD we 
readily obtain two possible sets for the edge parameters: 



There are three cases to consider. 

(i) Let u^ — u 2 , > 0. Then the second set (J90|) must be discarded by failure to satisfy the 
reality condition for the wavenumbers. The first set (fH9~|) satisfies both the reality condition 
and the "entropy" conditions. Thus, for u 2 + — u 2 _ > one gets the dispersive shock with 
the negative dispersion wave pattern. The width of the dispersive shock is (s + — s~)t = 
\{u 2 + — u 2 _)t. The amplitude of the lead soliton is found from the relationship between the 
velocity and the amplitude for the mKdV(d) solitons moving against the background u, 
which is obtained by an elementary analysis of the travelling wave solution in the soliton 
limit and has the form : c s = —u 2 ± 2ua s /3 — a 2 / 6. Different signs in this expression 
correspond to solitons of different polarity supported by the mKdV(d) equation. Now, after 
simple calculation one gets a + = 2(\u + \ — For u + u_ > this coincides with the KdV 




(89) 




(90) 
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result (|31|) (to make such a comparison in the case when u + > u~ > 0, the KdV equation 
should be taken with the negative sign for the nonlinear term, which would result in the 
change A — ► —A in (|31j)). The location of the mKdV(d) dispersive shock given by (|89p. 
however, is different from the respective KdV case as for the KdV equation with the same 
initial data one would have \s + — s~\ = §(|w+| — |w-|)£ rather than — u 2 _)t. 

(ii) u 2 + = u 2 _, which implies u + = w_ (trivial case) or u + = — w_. The width of the 
dispersive shock becomes zero, which implies that symmetric initial discontinuity with u + = 
— m_ does not break and propagates as a whole with the velocity — u\. This corresponds to 
the exact solution of the mKdV(d) equation in the form of a smooth kink (in the Whitham 
approximation its width is equal to zero, hence the propagating discontinuity). 

(iii) If u\ — u 2 _ < 0, the first set of parameters ()89|) fails to satisfy the reality condition for 
the wavenumbers while the second set (J9U|) does not pass the "entropy" test. This implies 
that there exists a genuine global solution to the dispersionless limit which does not require 
a dispersive shock. This solution is a rarefaction wave x/t = —u 2 confined to the interval 
— u 2 _t ^ x ^ — u 2 + t. 

Of course, these results could be inferred from the known periodic solution and the full 
modulation system for the mKdV(d) equation in Riemann invariants (see 42j for instance) 
but here they have been obtained without invoking the integrable structure mKdV(d) equa- 
tion, as an illustration of the effectiveness of the general transition conditions for a dispersive 
shock. 

5 Bi-directional dispersive hydrodynamics 

5.1 General setting and formulation of the problem 

We now generalise the obtained results to a physically important case of 2 x 2 strictly 
hyperbolic systems modified by weak dispersion. We represent such a system in a symbolic 
form 

4>t = K 4 (0, <j) x , 4> xx , 4> xt . . . ) . (91) 
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where <fi is 2-vector, K 4 is vector differential operator of fourth order with respect to spa- 
tial/mixed derivatives so that the system (f9~T]) has real linear dispersion relation u = u>o(k) 
so that ujq — > k as k <C 1 (weak dispersion). We assume that the system (JDT]) has at least 
three independent conservation laws 

d Pi dQi 

lf + 17 = - > = 1 * 3 - (92) 

For convenience of explanation we associate two conserving densities P12 with the "gas" 
density p and the momentum pu and assume that the dispersionless limit of (|91|) has the 
form of the gas-dynamic Euler equations for the isentropic ideal gas 

Pt + [pu)* = , (jm) t + (pu 2 + p(p)) x = 0, (93) 

where p(p) is the pressure in the corresponding gas dynamics. It should be emphasized that 
generally speaking the assumption about the gas dynamic 'core' of the system (}9~T]) is just 
a convenient (and in many cases physically relevant) way to convey our ideas - actually 
there is no need to restrict oneself with this particular form of the dispersionless limit. The 
only property of the dispersionless limit system which will be used below is that it can be 
represented in the Riemann form (which is always the case for quasilinear 2x2 systems). 
So actually one should put quotation marks for the "density" , "velocity" and "pressure" in 
this section. 

As earlier, we assume that system (J9~T|) supports the single-phase periodic travelling wave 
solutions for (f> = (p,u): 

(f)(x, t) = (f>(x - ct) , <p(9 + 2n/k) = <p(6) , (94) 

which are parametrised by four constants, say p, u, k, a, where p and u are the mean density 
and mean velocity respectively, and k and a are, as usual, the wavenumber and the amplitude. 
We assume that the potential curve G(X) in the travelling wave equation \ 2 e = —G(X) has 
generic properties outlined in Section 4 (see (|78jl. (JZHJ)) such that the solution (|9*l)l admits 
the limiting transitions to a linear wave as a — > and to a solitary wave as k —>■ 0. 

The described class of systems (19 ip is quite broad and includes some known integrable 
models such as defocusing nonlinear Schrodinger equation and Kaup-Boussinesq system 
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|44j . As physically important examples of bi-directional non-integrable systems that pos- 
sess the above general properties (including the gas dynamics form of the dispersionless 
limit) one can indicate the Green-Naghdi system for fully nonlinear shallow water gravity 
waves [13] and its multi-layer generalisations [21] , the systems for nonlinear ion-acoustic and 
magnetoacoustic waves in collisionless plasma (3B], [Tfi] . and many others. 

We consider initial data for the system (pTTj) in the form of a step for the variables p and 

u: 

t = : p = p~ , u = u~ for x < 0; p = p + , u = u + for x > , (95) 

where p ± and are constants. 

Analytical studies of the decay of an initial discontinuity problem in integrable dispersive 
wave equations (see for instance jUj, jHj, jS]) as well as direct numerical simulations for non- 
integrable systems ( see j2Ej j2H]- PB] and references therein) suggest that the asymptotic 
solution for the decay of an arbitrary initial discontinuity problem generally consists of 
three constant states separated by two expanding waves: centred rarefaction wave(s) and/or 
dispersive shock(s), which is quite natural taking into account the "two- wave" nature of the 
system 1)91)). The structure and the qualitative properties of the dispersive shock are the 
same as in the case of simple-wave lead dispersive equations (see Section 4) so we will model 
such a dispersive shock by the expansion fan solution of the corresponding Whitham system, 
which can be represented in the form 

d— d — 

—P j (p,u,k,a) + —Q j (p,u,k,a) = 0, j = 1,2,3, (96) 
d d 

+ —u(p,u,k,a) = 0, (97) 

where p and u are the density and velocity averaged over the family (|94)1. k is the wavenumber 
and a is the wave amplitude; also as follows from the formulation Pi = p, P 2 = ~pu. Using 
the arguments presented in Section 3.1, we postulate the following fundamental property of 
the Whitham equations (|9l)|) for weakly dispersive nonlinear wave equations ()91|) : 

The averaged "dispersive-hydrodynamic" conservation laws (J96|) admit exact 
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reductions to the dispersionless limit system (J93|) for a = or k = 0: 

Pt + {pu)x = , {pu)t + (pw 2 + p{p))x = 0, (98) 

so that the "energy" equation for P 3 becomes a consequence of (J98|) . Of course, in concrete 
instances this physically transparent property can be easily verified by direct calculation 
using the asymptotic behaviour of the potential curve G(X) for the nearly linear and nearly 
soliton configurations (see (fTHj). (|T9*j) ). 

The wave number conservation law (pTTj) in the linear limit assumes the form (cf. (^HJ)): 

a = 0, k t + (u (p,u,k)) x = (99) 

where = co>o(p, u, k) is the linear dispersion relation for the system (|91j) obtained by consid- 
ering an infinitesimally small perturbation of the ground state p = p, u = u. Owing to the 
two-wave nature of the system (}9~T]) there are two branches of the linear dispersion relation. 
To avoid unnecessary complication, it will be assumed that, unless otherwise specified, we 
consider the branch corresponding to the right-propagating waves. 

Now we can formulate the natural boundary conditions of the Gurevich-Pitaevskii type 
for a bi-directional Whitham system ()9(i)l . (|9Tj) . Again, for simplicity of presentation we 
assume the negative dispersion wave pattern and later will remove this restriction in the 
formulation of the transition conditions. In the decay of an initial discontinuity problem 
we are interested in the similarity solutions of the modulation equations so we present the 
natural matching problem for modulation variables at once in the form analogous to (|3l)|) 

x = s~t : a = 0, p = p~ , u = u~ , 

(100) 

x = s + t : k = , p = p + , u = u + , 

where the dependencies of the edge speeds s ± on the initial jump parameters p + , p~ ,u + , u~ 
are to be found. 

It is clear that the similarity solution of the fourth order Whitham system is parametrised 
by only three constants. Therefore, one must impose an additional restriction on the values 
of the parameters p ± , to select the family of admissible jumps across the dispersive shock. 
Of course, the necessity of an additional jump condition for the dispersive shock generated 
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in the decay of an initial discontinuity is inherently implied by the two-wave nature of the 
system under consideration. Finding this restriction is equivalent to the extracting the family 
of the initial discontinuities resolving into a single dispersive shock propagating in a given 
direction (we associate the direction of the dispersive shock propagation with the direction 
of the corresponding linear characteristic). In dissipative gas dynamics, such admissible 
discontinuities are selected by the Rankine-Hugoniot conditions following from the balance 
of mass, momentum and total energy across the shock. So, before we proceed with the 
determination of the dispersive shock edges as we did in the simple-wave led equations, our 
task is to obtain a replacement for the classical shock curve in the bi-directional dispersive 
hydrodynamics which would have a form of a relationship 



We note that in the conservative dispersive shock transition, the pressure p is the function of 
the density only (the thermodynamic entropy does not change in the course of the dispersive 
shock propagation) so, in contrast to the classical shock theory ^Hj, the pressure should not 
appear in the transition relation (jlOlj) . 

5.2 Dispersive shock curve and "entropy" conditions 

We proceed with the dispersive shock curve using a simple method recently proposed in jlH] • 
The idea is, using the conservative nature of the dissipationless dispersive shocks to "reverse" 
the dispersive shock solution in time. It is then "natural" to expect that the dispersive shock 
would convert into a rarefaction wave of the dispersionless Euler system (193(1 . One can then 
"extract" the required transition relation from simple analytic solution for this rarefaction 
wave which is easily constructed using Riemann invariants (see for instance [HI], P4J). Of 
course, such an intuitive argument requires mathematical justification since it is not obvious 
why the "inverse" resolution should happen through a single rarefaction wave. 

Our construction is based on the following geometric consideration in the characteristic 
(x, t)-plane. First we represent the Euler system (|D3J) in the Riemann invariant form 



$(p + ,p ,u + ,u ) = 



(101) 



r t + V r (r,l)r x = 0, 



k + V l {r,l)l x = 0, 



(102) 
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where 



r = u + 




= u — 




^V; V r ,i = u±<r(p). (103) 



Here p — const and er(p) = {dp/ dp) 1 / 2 is the "sound speed". The initial conditions (|93j) are 
rewritten in terms of the Riemann invariants as 



where r^ 1 = r(p ± ,w ± ), I = l(p ,u ). The corresponding four characteristic directions at 



We consider the "right-propagating" dispersive shock, which occurs when the characteris- 
tics of the r-family intersect i.e. when V~ > . Then, assuming the modulation description 
of the dispersive shock with the aid of the expansion fan solution of the Whitham equations, 
one can observe that, due to (assumed) hyperbolicity the solution of the GP problem (j9*Bj) . 
(|9*7|). ()100|) can be, in principle, constructed geometrically using characteristics. The continu- 
ity matching (jlOOj) for the mean values in such a construction is replaced with the equivalent 
continuity matching for the Whitham and Euler characteristics (see jHH])- One can also 
observe that although such a construction can be very complicated for t > 0, it can be 
readily realised for negative t. Indeed, since the Whitham expansion fan has zero width at 
t = (hyperbolicity) the initial conditions ()104|) are specified entirely in the external, "dis- 
persionless" domain of the space-time of the GP problem. Thus the global solution of the 
GP problem can be continued backwards along the characteristics of the 2x2 dispersionless 
limit equations (jl(J2|) . By construction, this solution must (i) be three-parametric (see (jlUlJl ) 
and (ii) satisfy the inequality V~ > at t — 0. 

There is a unique three-parametric solution to the hyperbolic system ()102)) in the lower 
(xt) half-plane satisfying the described restrictions. This solution is a centred expansion fan 
(in (—x, —t) - coordinates) given by the expressions 



r = 



r ,1 = 1 for x < ; 



r = r , 1 = 1 for x > 



(104) 



t = are V± = K-(r ± , / ± ) and Vf = V^r* / ± ). 



I = l Q = constant 



(105) 
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r = r~; x > a 2 t; 

V r (r,l ) = x/t; < x < a 2 t; (106) 

r = r + ; x < ait. 

Here 

ai = V r + , a 2 = V7~ , a 2 > ai . (107) 

The solution (|105J) - ()107j) is characterised by three parameters (r + ,r~,Zo) and if exists for 
all t < 0. Since this solution represents a continuation, along the characteristics of the full 
solution to the GP problem, one can regard this solution evaluated at a fixed moment to < 
as new initial conditions for the same GP problem. Therefore, one can extract from it the 
global restriction (jlOlj) imposed on possible values of the hydrodynamic variables at the 
opposite sides of the dispersive shock. This restriction follows from f)105j) and has the form 

r = l + , (108) 

which in view of the relationships (jK)Hj) can be represented in an explicit gas-dynamic form, 

u- -u + = f ^-dp. (109) 

Jp+ p 

Given the state in front of the dispersive shock (p + ,u + ) this relation yields all admissible 
states (p~,u~) behind it i.e. it represents the equation of the p-u diagram of the dispersive 
shock. One can see that the whole above construction is subject to the additional inequalities 



Vf < s~ < V~ , V r + <s + , s + > s- . (110) 

These inequalities ensure fulfillment of our original requirement about the single-wave res- 
olution of a step (i.e. three-parametric solution). Indeed, the number of parameters char- 
acterising the solution of the hyperbolic system in some domain is equal to the number of 
families of characteristics transferring given initial or boundary data into this domain (see 
for instance jlH])- The inequalities (jllOj) require that only three of the four families of clas- 
sical characteristics x/t = (namely, x/t = V~ , x/t = V T + , and x/t = V^) transfer 
initial data ()104|) from the x-axis into the dispersive shock domain (see Fig. 4). One can see 
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that the inequalities ()110|) represent an extension of the "entropy" conditions (f75|) to the 
bi-directional case. 

a) 




x='a l t 




Figure 4: Qualitative behaviour of characteristics in the "simple-wave" decay in dispersive 
hydrodynamics. Broken lines: the dispersive shock boundaries. Dotted lines: the "mirror" 
expansion fan boundaries, (a) Families dx/dt = transfers values of r, (b) Families 
dx/dt = transfers values of /. 

Now we discuss briefly the meaning of the obtained transition relation (J109J) . One can 
see that it coincides with the relationship between any two pairs p, u in the simple wave 
solution of the isentropic gas dynamics (see for instance [49J). A nontrivial fact is that the 
corresponding wave is the simple wave of compression which breaks after a certain time 
interval. Contrastingly, the similarity solution of the Whitham equations satisfying the rela- 
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tionship ()105|) and describing the expanding dispersive shock does not break. This solution 
essentially represents a compression fan, which does not exist in classical gas dynamics. It 
is natural to call the dispersive shock satisfying the relationship (jl05|) a simple dispersive 
shock. We emphasize that, according to the matching conditions (jlOOj) the relationship (jl09|) 
is only valid for the boundary values of the modulation parameters p and u and, of course, 
does not hold within the dispersive shock region. 

In the solutions of the GP problem for integrable systems, the simple dispersive shock 
condition is a mere consequence of the constancy of one of the Riemann invariants of the 
Whitham system (see for instance). In non-integrable case, when the Riemann 

invariants are not available for the Whitham system, the condition f)109j) is not obvious at 
all. Also, one should keep in mind that inequalities 1)110)1 represent a necessary part of the 
transition conditions. 

We recall that the above consideration was concerned with the dispersive shocks prop- 
agating to the right. For the left-propagating dispersive shocks an analogous system of 
relations would include the zero jump condition for the classical Riemann invariant r in- 
stead of I in (1105)) and the "entropy" inequalities analogous to (jllO)) would have the form 

V+ <s + < V r + , s~ < Vf , s + > s" . (Ill) 

The simple dispersive shock curve in the form (J109j) has been proposed for the first time, 
using intuitive physical arguments, by Gurevich and Meshcherkin in the context of 
collisionless shocks in plasma. Later, this condition has been interpreted in [33J in terms of 
the "local Riemann invariant" transport through the Whitham zone. 

We note in conclusion that, in a more general case, when the dispersionless limit has the 
form other than isentropic gas dynamics, one should use the simple dispersive shock curve 
in the invariant form ()108j) . 
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5.3 Speeds of the dispersive shock edges 

We now proceed with the dispersive shock edges similarly to Sections 3.2, 3.3. We consider 
a characteristic equation for the Whitham system (jHSJ), (j97j) 

bfdip + bfdiU + bfdik + bfdta = , (112) 

where b^(p, u, k, a) is the left eigenvector of the coefficient matrix A(p, u, k, a) of the modula- 
tion system (|96|h (|97p and i is a number of the characteristic. Let a = on this characteristic. 
Then for the remaining three variables we have the ordinary differential equation 

bf (jo, u, k, 0)dip + bf (p, u, k, 0)diU + b^ (p, u, k, 0)dik = (113) 

which necessarily is a characteristic equation for the reduced as a = modulation system 

(H3), dSZD: 

a = 0, p t + (pu) x = 0, (pu) t + (pu 2 +p(p)) x = 0, k t + {u {p,u, k)) x = 0. (114) 

Now we impose an additional constraint 

F(p,u) = Co, (115) 

where Cq is a constant. It will be shown later that this constraint leads to a self-consistent 
set of the transition conditions. One can see that Eq. ()115j) is consistent with the reduction 
(|114j) and, therefore, is compatible with (|113|) iff the function F(p,u) is one of the Riemann 
invariants l,r defined by ()103|) . Setting Eq. ()115|) into Eq. ()113p we obtain a 2-form which 
can be integrated to give the local relationship between admissible values of p and k on the 
characteristic, 

Mu,k;C ) = d, (116) 

C\ being a constant of integration. 

Alternatively, the relationship ()116|) can be obtained by a direct substitution of Eq. ()115j) 
into the system ()114)) and then, by looking for the solution in the form k(u). Substitution 
of Eq. (fUBJ) into JHl yields (cf. flUJ)): 

a = 0, u t + V(u)u x = 0, kt + (u (u, k)) x = , (117) 
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where (see Eq. (1103(0 

V(u) = V r (p(u),u) , uj {u,k) = u (p(u),u,k) . (118) 

Thus, the problem essentially reduces to the simple-wave led case considered in Sections 3, 4 
and, therefore, yields the ordinary differential equation for k(u). In our case, however, 
this equation is additionally parametrised by the unidentified (yet) function F(p,u) and a 
constant Co in (jllf>|) defining the dependence p{u) in ()118|) . These can be found by applying 
the matching conditions (jl()()|) considered in the four- dimensional space of the field variables 
p, u, k, a. Since the integral ()115|) does not contain a and k it must apply to both boundaries 
a = and k = 0, which requires F(p~,u~) = F(p + ,u + ). Then, to be consistent with the 
dispersive shock curve (jl()8|) one must identify F(p,u) = l(p,u), which immediately yields 
Co = l~ = l + ■ Now, the trailing edge speed s~ is found with the aid of Eq. ()84a ). 

The leading edge edge is handled in exactly the same way, by reducing the "hydrody- 
namic" part of the modulation system for k = to the simple wave equation with the aid 
of the constraint ()115|) . Then one proceeds with the conjugate variables as in Section 3.3 
to get a general simple-wave analog of the ordinary differential equation (joTlj) i.e. Eq. (j%3l 
for k(u). At last, using the dispersive shock curve (jl()8|) the function F(p,u) in (J118)) is 
identified with the same classical Riemann invariant /, which justifies self-consistency of the 
whole construction. Then the leading edge speed is found via the conjugate wavenumber k 
by Eqs. (JMb ). (|H3~|) where all the necessary ingredients are given by Eq. ()118|) . 

The whole construction is subject to the "entropy conditions" ()85|) which should be 
complemented by an additional inequality s~ > Vf (see ()110J0 occurring due to the presence 
of the fourth characteristic family in the bi-directional problems (see [IE] for an example 
where this inequality is violated). 

It is clear that analogous formulas can be obtained for the left- propagating dispersive 
shock: one should just replace in Eq. ()118|) V r (u) with Vi{u) and /(p, u) = C with r(p, u) = 
Cq where Cq = r~ = r + according to the corresponding dispersive shock curve. Also, one 
should use the "left" branch of the linear dispersion relation and the "entropy" conditions 

(EU). 
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Some additional details of the derivation of the trailing and the leading edge curves for 
bi-directional systems can be found in |48| . 

5.4 Example: fully nonlinear ion-acoustic collisionless shock in 
plasma 

As an example of an effective construction of the dispersive shock conditions in a non- 
integrable system we consider classical system of equations describing finite-amplitude ion- 
acoustic waves in a two-temperature (T e ^> Tj) collisionless plasma (see for instance ^Hj) 

Pt + {pu) x = 0, 

u t + uu x + if x = 0, (119) 
Vxx = e v - p . 

Here p and u are the ion density and velocity and (p is the electric potential; all dependent 
variables are dimensionless. We note that the system does not contain the time derivative 
of ip so p is not a "real" dependent variable as concerns the 2x2 representation 1)91)1 . A 
direct numerical simulation of the decay of an initial discontinuity in Eqs. ()119)) has been 
performed in f2Jo\ . 

The system (1119)1 supports periodic travelling waves allowing for linear and solitary wave 
limits and also possesses (at least) four conservation laws [38J. Thus we can apply the 
methods developed in this paper for obtaining the dispersive shock conditions. We consider 
the right-propagating simple dispersive shock. 

In the dispersionless limit tp = In p, which yields the Euler isentropic gas-dynamic equa- 
tions ()93p with the equation of state p(p) = p. The corresponding Riemann invariants and 
characteristic velocities are (see 1)102)1 . I) 103)1 ): 

l — u — lnp, r = u + \np, V\ — u — 1 , V r — u + l, (120) 

while the linear dispersion relation for the right-propagating modulated waves has the form 

u (k, p, u) = k[u + (1 + k 2 /p)-^ 2 ] . (121) 
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The dispersive shock transition curve (jl09j) then assumes the form u~ — u + = ln(p~/p + ). 
Without loss of generality we put u + = 0, p + = 1. Then the relationship between p and u 
in Eq. (jll8j) becomes u = In p. As a result, we get all the necessary ingredients for the basic 
ordinary differential equations 



V{u)=u + 1, to (u,k) = k[u + (l + k 2 /e a )- 1/2 }. (122) 



Then the equation ()82|) after elementary transformations, assumes the form with sepa- 
rated variables 

d ^ = - , (1 + a) \ , a(l) = l, (123) 
du 2(1 + a + a 2 ) ' K ' K J 

where a = (1 + k 2 je u )~ l l 2 . Integrating ()123|) we get 

M + 21na + ^^ = 0. (124) 
1 + a y ' 

Now, using the formula (JS4*k ) we obtain a simple equation determining velocity of the trailing 
edge s~ in terms of the density ratio across the dispersive shock d = p~ / p + = p~ , 

2 (q- — In d) 1 / 3, — 1 

tad+ - H s-- tad) J J —J jw - i . ( 125 ) 

The leading edge is handled in a completely analogous way. The solitary wave dispersion 
relation is obtained from the linear dispersion relation ()122j) as uj s (u, k) = —iuo(u,ik) which 
is to say 

u s {u, k) = k[u + (1 - k 2 /e a )- 1/2 } . (126) 

Then, integrating ()83|) we obtain k{u) (it is also convenient to introduce a = (1 — k 2 /e u )~ l l 2 
as an intermediate variable instead of k). Setting it into (|84b) we eventually get the equation 
for the leading edge 

s+ - 1 

21ns + -— = \nd. (127) 

It can be easily verified that the "entropy" conditions ()110j) are satisfied for all values of 
the density ratio across the dispersive shock. The curves s + (d) and s~(d) are presented in 
Fig. 5 and demonstrate very good agreement with the results of direct numerical simulation 
of the decay of an initial discontinuity for the system ()119|) presented in (see Fig. 6 in 
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Figure 5: Speeds of the edges of the ion- acoustic dispersive shock transition versus density 
ratio d — p~ j p + across the shock; upper graph: s + - leading edge, lower graph: s~ - trailing 
edge 

this paper analogous to our Fig.5). From theoretical point of view such an agreement can 
be regarded as a strong indication of validity of the modulation theory in a certain class of 
non-integrable initial value problems where rigorous derivation of the Whitham asymptotic 
as a zero-dispersion limit is not available. 

The weakly nonlinear asymptotic decompositions of (jl25J) and (jl27J) for small 5 = d— 1 <C 
1, have the form 



matter of fact, agree with the Gurevich-Pitaevskii solution for the KdV solution (|3U|) 
which is another confirmation of validity of the obtained modulation solution. Curiously, 
as is clearly seen from Fig. 5 , the fully nonlinear dynamics of the leading (solitary wave) 
edge s + is very well approximated by the weakly nonlinear asymptotic (|128|) in a broad 
range of the density ratios d while the speed of the trailing (harmonic) edge s~ demonstrates 
significant qualitative and quantitative deviations from its weakly nonlinear counterpart even 
for quite moderate values of d. 

A detailed comparison of the modulation transition conditions obtained here with results 



s 
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of direct numerical simulation for the (non-integrable) Green - Naghdi system describing 
fully nonlinear shallow- water waves will be published in [5T] . 

6 Some restrictions 

We first address the accuracy of the obtained results in the context of the original (disper- 
sive) system and outline some restrictions of our analysis of the dispersive shocks. Our main 
assumption was about modelling the a dispersive shock with the aid of the expansion fan 
solution of the Whitham equations. The Whitham equations, being obtained by an asymp- 
totic procedure, have inherent accuracy restrictions on their applicability to certain wave 
regimes. In our (GP) formulation, in contrast to the formal zero-dispersion limit approach, 
the small dispersion parameter determining the accuracy of the Whitham approximation is 
assumed to occur in the solutions and is defined by the ratio of the typical wavelength to 
the characteristic scale of the modulation variations (in our case, the width of the dispersive 
shock L ~ t), i.e. 

e~(Jfet) _1 «l. (129) 

which implies that that modulation description of the dispersive shock is valid asymptotically 
as t ^> 1. The Whitham method of averaging conservation laws jTHj also requires that the 
relative variations of modulation parameters over the characteristic wavelength scale ~ k~ l 
be 0(e) which yields the criterion 

\k x \k' x ^r l . (130) 

One can now observe that, in the modulation solution for the dispersive shock, for any t, 
however large, there is a certain vicinity of the of the leading edge where both criteria ()129|) . 
(|13Uj) are violated due to small (and rapidly changing) values of k. This is basically a re- 
flection of non-uniformity occurring as k — > in the formal perturbation decompositions 
(see e.g. [20]) equivalent to leading order to the Whitham method. So a natural question of 
applicability of the modulation description to the dynamics of the leading edge of the dis- 
persive shock arises. While this question does not appear in the the rigorous zero-dispersion 
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limit approach |TTj in which the Whitham equations are derived as a certain asymptotics 
in the initial-value problem and e does not depend on t, it has to be addressed in our case 
when the Whitham description is assumed on heuristic foundations in the absence of rich 
integrable structure. 

We first determine the behaviour of the function k(x) in the vicinity of the leading edge of 
the simple dispersive shock for fixed t. It is instructive to make first an estimate for the KdV 
case, which can be made by using the wave number conservation law in the self-similar form 
—sk' +u' — 0, s — x/t. Now setting the asymptotic decompositions for k and w as m' < 1, 
which are readily obtained from @ we get for the vicinity of the leading edge s + -s C 1 the 
following asymptotic behaviour: k ~ fco/hi^" 1 , u—ks + ~ k5, 5 ~ {s + — s)(ln l/(s + — s)) -1 . It 
is clear that this behaviour is automatically generalised to other nonlinear dispersive systems 
for which the potential curve G(u) in the travelling waves has the asymptotic behaviour (J7%|) . 
(|79l) in the nearly linear and nearly soliton configurations. 

Now, using the obtained asymptotics for k, we estimate the relative widths a = Ax/L of 
the vicinities of the leading edge x + = s + t where the criteria ()129j) . (jlHOjl are violated. For 
(THZ91 we get a x ~ exp(-i) while for (fHUl) we get a 2 ~ t" 1 hit, i.e a 2 > <J\. Still, a 2 (t) -> 
as t — > oo and, therefore, the Whitham description of the dispersive shock is asymptotically 
as t ^> 1 valid for all s~t < x < s + t. 

The modulation approach used in this paper requires the existence of the single-phase 
periodic solutions. For non-integrable systems such solutions often exist only within a certain 
domain of parameters. Typically the role of the "critical" parameter is played by the wave 
amplitude, so that for a > a cr the periodic solution (or a solitary wave) ceases to exist 
because of occurrence of breaking or cusp- type singularities. This usually violates the "single- 
flow" type assumptions used in the derivation of the original system. For instance, for 
ion-acoustic system (|119|) the critical solitary wave amplitude for the potential tp is a cr w 
1.3 For a s > a cr the "hydrodynamic" single-flow description of the two-temperature 
plasma becomes not applicable and a more general, kinetic plasma theory should be used. 
The critical value for the density ratio across the ion-acoustic dispersive shock A cr can 
be found by setting the critical value of the solitary wave speed c s (a cr ) ~ 1.6 J2S1 into 



49 



Eq. (|127j) : s + (A cr ) pa 1.6, which yields A cr pa 2.0. Another possible restriction on the 
application of the obtained transition conditions partially addressed in Section 4 is imposed 
by the requirement of the modulational stability of the dispersive shock solution. For some 
equations the type of the modulation system can change from hyperbolic to elliptic depending 
on the initial data. One of the examples is provided by the so-called Kaup-Boussinesq 
(integrable) system describing bi-directional shallow water waves where the linear dispersion 
relation yields complex frequencies if the wavenumber exceeds a certain value |44j . The 
corresponding restrictions on the jump values across the dispersive shock can be found from 
the conditions (jHEJ)- 

7 Conclusions 

A new method has been proposed to analyse the dispersive shock transition for a broad 
class of weakly dispersive nonlinear wave equations. No assumptions of integrability have 
been made which, in particular, allows one to apply the developed method to problems of 
fully nonlinear dispersive wave dynamics. The transition conditions have been derived by 
assuming the description of the dispersive shock with the aid of the expansion fan solutions 
of the associated modulation (Whitham) equations. The Whitham system was assumed to 
be hyperbolic for the solutions of our interest, which seems a reasonable assumption for the 
outlined class dispersive systems (at least for a certain type of initial data). This assumption 
can also be viewed as an inference from direct numerical simulations. 

The analysis has been performed using the generalisation of the so-called Gurevich- 
Pitaevskii problem formulated originally for the averaged KdV equation. It has been shown 
that the Gurevich-Pitaevskii type natural boundary conditions for the mean flow in the dis- 
persive shock region can be translated into information about degeneration of the families of 
characteristics along the dispersive shock boundaries. The latter implies certain restrictions 
on admissible values of the hydrodynamic and wave parameters at these boundaries. It has 
been shown that, for the problem of an initial step resolution, these restrictions can be found 
without presence of the Riemann invariant structure for the Whitham system. 
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As a result, we have derived a complete set of transition conditions for dispersive shocks 
linking two different constant hydrodynamic flows. The analysis has been performed as for 
single-wave equations so for bi-directional systems. The transition conditions have been 
derived in a general form and include: i) a "dispersive" analog of the traditional jump 
conditions; ii) the equations for the boundary curves of the dispersive shock; hi) a set of 
inequalities similar to the entropy conditions of traditional gas dynamics. 

Remarkably, the whole set of the transition conditions is constructed in terms of the 
linear dispersion relation and the dispersionless nonlinear characteristic velocities of the 
system under study. The obtained conditions allow one to "fit" a dispersive shock into 
the classical dispersionless solution without complicated analysis of the internal structure 
of the dispersive shock similarly to classical shock theory in traditional gas dynamics. The 
method can be useful not only for the analysis of non-integrable dispersive systems where full 
solution is not available but also for integrable systems when one is interested only in main 
physical parameters of the dispersive shock transition. As an additional bonus, the developed 
method allows for obtaining the lead solitary wave amplitude, the major parameter observed 
in experiments. 
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